Magnetic field generation from cosmological perturbations 
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In this letter, we discuss generation of magnetic field from cosmological perturbations. We consider 
the evolution of three component plasma (electron, proton and photon) evaluating the collision term 
between elecrons and photons up to the second order. The collision term is shown to induce electric 
current, which then generate magnetic field. There are three contributions, two of which can be 
evaluated from the first-order quantities, while the other one is fluid vorticity which is purely second 
order. We estimate the magnitudes of the former contributions and shows that the amplitude of the 
produced magnetic field is about ~ 10 _19 G at lOMpc comoving scale at the decoupling. Compared 
to astrophysical and inflationary mechanisms for seed-field generation, our study suffers from much 
less ambiguities concerning unknown physics and/or processes. 
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INTRODUCTION 

There are convincing evidences that imply existence 
of substantial magnetic fields in various astronomical ob- 
jects. Not only galaxies but systems with even larger 
scales, such as cluster of galaxies and extra-cluster fields, 
have their own magnetic fields (for a review on cosmo- 
logical magnetic fields, see e.g., |l|). Conventionally the 
magnetic fields in galaxies, and possibly in clusters of 
galaxies, are considered to have been amplified and main- 
tained by dynamo mechanism. However, the dynamo 
mechanism needs the seed magnetic field and does not 
explain the origin of the magnetic fields. 

There have been many attempts to generate the seed 
fields. One of the approaches to this problem is to gen- 
erate magnetic fields astrophysically, often involving the 
Biermann mechanism This mechanism has been ap- 
plied to various systems: large-scale structure formation 
jjj, ionizing front protogalaxies [j| and supernova 
remnant of the first stars |(|. These studies showed the 
possibilities of magnetic-field generation with amplitudes 
of order 10 -16 ~ 10~ 21 G, which would be enough for the 
required field ~ 10~ 20 G. 

On the other hand, cosmological origins, which are of- 
ten concerned with inflation, can produce magnetic field 
with coherence lengths of much larger scales, which is 
typically the horizon scale or possibly super-horizon scale 
S E3 03 [HI Also they can often produce fields with 
a wide range of length scales, while astrophysical mech- 
anisms can produce fields only with their characteristic 
scales. For a constraints on magnetic field with cosmo- 
logical scale, see [l2|. 

In this Letter, we consider magnetic-field generation 



from cosmological perturbations after inflation. Around 
and after the decoupling, coupling among charged par- 
ticles and photons become so weak that electric current 
can be induced by the difference in motions of protons 
and electrons. This electric current leads to generation of 
magnetic fields. It is well known that vorticity of plasma 
produces such electric current [lH HI HI uR Il7l [Tsj . 
However, because vorticity is not produced at the first 
order in cosmological perturbations, we must study the 
second order. We will consider equations of motion for 
protons, electrons and photons separately up to the sec- 
ond order, although equation of motion for photons does 
not appear explicitly. To study electric current appro- 
priately, we must treat the three components separately 
|l7| . Furthermore, we will evaluate the collision term 
between electrons and photons, which is dominant and 
essential for the magnetic-field generation, up to the sec- 
ond order. From the equations of motion for protons 
and electrons, we obtain a generalized Ohm's law, which, 
combined with the Maxwell equations, leads to an evo- 
lution equation for magnetic field. Our study is based 
on the cosmological perturbation theory, which is highly 
successful in the anisotropies of the cosmic microwave 
background, and suffers from much less ambiguities con- 
cerning unknown physics and/or processes compared to 
astrophysical and inflationary mechanisms for seed-field 
generation. Thus our results are robust. 
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FORMULATION 

Euler equations for proton fluid and electron fluid are 
given by 

(s\ + «V) (r>% + r^) = eg + qg> (i) 

(S\ + (t^ + T^j) = + C^>\ (2) 

where T^i-i and T-£ Mp ,>. are the energy-momentum tensor 
of proton (electron) fluid and electromagnetic field cou- 
pling to protons (electrons) current, respectively. Here 
/i, v = 0, f , 2, 3 and i = 1,2, 3. The projection of the di- 
vergence of the energy-momentum tensors are computed 
as 

= (p + p)u»u\^ + (6\ + u l u,)P^, (3) 
(S^ + u'u^T^^-fF^, (4) 

where p, P and are the energy density, pressure and 
electric current, respectively. The r.h.s. of Eq. 0) and 
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(J2J represent the collision terms. C pe = —C ep is the 
collision term for the Coulomb scattering between pro- 
tons and electrons. This term leads to the diffusion of 
magnetic field and can be neglected in the highly con- 
ducting medium in early universe. On the other hand, 
the collsion terms for the Thomson scattering between 
protons (electrons) and photons are expressed as cj|i. 
Because the collision term for the protons can be ne- 
glected compared to that for the electrons, difference in 
velocities of protons and electrons will be induced which 
leads to electric current. This electric current becomes a 
source for magnetic field. 

Now we evaluate the collision term for the Thomson 
scattering: 

7(P<) + -» l{Pi) + e"(«0, (5) 

where the quantities in the parentheses denote the par- 
ticle momenta. The collision integral Cye [/(Pi)] f° r this 
scattering is given as 



2tt 4 r d 3 p' r d 3 q r d 3 q ' 

~p~ J (2tt) 3 2E 1 ( P / ) J (2n) 3 2E e (q) J {2n) 3 2E e {q') 

x \M\ 2 6 [E, ( (p) + E e {q) - E^p') - E e (q')} S (3) [ Pl + ft - v\ - ft'] {/ 7 - / 7 (R>/e(ft)} 

n 1 ^\Mnm)-Mp t )}\s { p-^ + (p t - P> : d J^ 



Am*p J (2ir) 3 p 



dp 1 



(6) 
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where -E 7 ( e ) and / 7 ( e ) are the energy and the distribution 
function of photons (electrons), respectively, \M\ 2 is the 
scattering amplitude and m e is the electron mass. Then 
we obtain the collision term in the Euler equation (0) , as 



|^Cg>[/(pO] 



A7) 



where <jt is the cross section of the Thomson scattering. 
Here moments of the distribution functions are given by 



d 3 p 
(2tt) 3 

J^P = ^P 7 U 7 ' 

d 3 p 



Pf-yiPi) = Pi, 

4 
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3 P l fe(Pi) = PeU\ 
d 3 p _ n , , . , , 1 „;; I 



(27T) 



(8) 
(9) 
(10) 



( . >7r)3 -p-W/7(Pi) = qP^J + §P7^> (11) 



where is photon anisotropic stress. It should be 
noted that the collision term (0 was obtained non- 
perturbatively with respect to the cosmological pertur- 
bation. 

Altogether, the Euler equations for protons and elec- 
trons are written as 



enu^F; = 0, 



AaTP^n 



(12) 
(13) 



where m p is the proton mass, and the pressure of proton 
and electron fluids are neglected. We also assumed the 
charge neutrality: n — n e ~ n p . Subtracting Eq. I|12|) 
multiplied by m e from Eq. I|13|l multiplied by m p , we 
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obtain 



m p m e 



J' 



+en(m p + m e )u^F 1 ^ - (m p - m e )j fl F i fi 



Am p <jT p-yTi 



1 



(14) 



where u M and j M are the center-of-mass 4-velocity of the 
proton and electron fluids and the net electric current, 
respectively, defined as 



f = en(u p l 



- m e 
<)■ 



(15) 
(16) 



where e yfc is the Levi-Civita tensor and B % = e v Fjk/2 is 
comoving magnetic field. We can expand the photon en- 
ergy density, fluid velosities and photon anisotropic stress 
as 



(0) 

Pi =P~r 



(i) 
P~ 



u = 1- 



(2) 

u 



(1) 

~-Ui 



(2) 
U i 



■ ■ , ^ =n« + • • • , (20) 

where the superscripts (i) denote the order of expansion. 
Remembering that B l is a small quantity, we see that 
most of the terms in Eq. I|19|) . other than the first and 
third terms, can be neglected. Thus we obtain 



Noting the Maxwell equations F£" — j M , we see that 
the first term in the l.h.s. of Eq. (|14fl are suppressed, 
compared to the second term, by a factor 



3 x 10 



-47 



/IMpc 



L 



(17) 



where c is the speed of light, L is a characteristic length 
of the system and uj p = y/ '4irne 2 /m e is the plasma fre- 
quency. Second order vector perturbations are contained 
in the covariant derivative of the electric current and were 
evaluated in [J^. They obtained the current magnetic 
field of order 10 _29 G at IMpc scale. As we will see be- 
low, this is much smaller than a value we obtain in this 
paper, which justifies the neglection of the first term in 
the l.h.s. of Eq. fity. 

The third term in the l.h.s. of Eq. ifT^) is the Hall 
term which can also be neglected because the Coulomb 
coupling between protons and electrons is so tight that 
\u l \ ^> \u l — u\\. Then we have a generalized Ohm's law: 



B' 



c yk c 



3,k 
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4cr T P 1 Ajk 
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(i) 

Pi,k /<i) 
— — u 

(0) 
Pi 



/(l) (1) \ /(2) (2) \ 



1 f (l) (1), (l) (1), 

ueiM n l 7j + u el rr 



lj,k 



(21) 



where the dot denotes a derivative with respect to the 
cosmic time, and we used the fact that there is no vortic- 
ity in the linear order: e y& u < & = 0. The contributions 
of the first two terms in Eq. (|21|l were first noticed in • 
From this expression, we see that magnetic field cannot 
be generated in the linear order. Here it should be noted 
that the velocity of electron fluid can be approximated 

to the center-of-mass velocity at this order, u\ ~tt\ 

The linear-order quantities have only scalar compo- 
nents so that we can write as 



3e 



C\ (18) 



Now we derive the evolution equation for magnetic 
field. This can be obtained from the Bianchi identities 

F[»v,\] = 0, as 



K 

n 



J U 



-(k l kj 



(22) 
(23) 



= e ij '*W bM 



= vF-B 1 



n £ °J U ,k + £ L/ .7,fe 



Also we define p 7 



(0) 

= P 



and 5~, 



(i) 

= P 7 



IPi 



Then we 



can rewrite the Fourier component of B l in terms of the 



-{u % jB J — v? -B l ) H ^-(B 3 u l — B l u J ), (19) Fourier components of <5 7 , u, u-y and n 7 , as 
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B\K,t) = I <II',,,U') 



d 3 k' 



k'„K k 



(2tt) 3 1 

^ J'dt'p^t') [^(K^-n^KX)] 



- %K ■ h 1 — 9k' 2 
8 7 {\K -k'\,t')S Ubj (k',t') - ^ -, u(\K - k'\, t')n 7 (fc', t') 



2Ak' 2 \K-k' 



(24) 
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FIG. 1: Power spectra (multiplied by k) in the source term of 
Eg. 1241 at z ~ 1100. The solid and dashed lines correspond 
to power spectra of (|<5 7 (« — it 7 )|) and (|n 7 w|), respectively. 
The spectra have been normalized to agree on large scales at 
present epoch. 

where ki is the unit vector in direction of ki, Su = u — 
Uj. Vorticities, which come from the second-order vector 

perturbation, are defined as fl l — e yfc SI k,j and fi* = 

e yfc ( u~ ( k.j- Eq. ((52} , which describes the evolution of 
magnetic field, is one of our main results. 

EVALUATION OF SOURCE TERMS 

Here we briefly evaluate the first two terms in the r.h.s. 
of Eq. (|24|l which are constructed from the first-order 
quantities. The evaluation can be easily accomplished 
by using linearized Boltzmann code such as CMBFAST 
|2fJ . In figure^we show the power spectra of these terms 
multiplied by k, kP$ u (k) and kPn, u (k), around the last 
scattering epoch z 1100. Here the power spectra are 
defined as 

(\5 7 (k){u- Uj }(k')\) =P Su (k)6W(k-k') (25) 
(|n 7 (fc) U (fc')|) = Pn u (fc)^ 3 )(fc-fc')- (26) 

For this illustration all cosmological parameters were 
fixed to the ACDM values, i.e., (h, n s , Qbh 2 , £ljs_, A) = 
(0.71,0.99,0.022,0.74,0.83). 

We found that these two terms contribute to the gen- 
eration of magnetic field at almost the same order of 



magnitude at all scales. This can be understood since 
<5 7 ~ it 7 (acoustic oscillations), u — u 7 ~ eRu, and 
II 7 ~ euj in the tight coupling approximation of baryon- 
photon fluid, where e being tight coupling parameter and 
R = 3pb/4py ~ C(l) around the last scattering epoch. 
From the spectra we can estimate the field strength 
of generated magnetic field at the decoupling epoch as 
B ~ ?^^k 3 P Su (k) - 10" 19 G at lOMpc comoving 
scale. If the magnetic field decays adiabatically accord- 
ing to the cosmic expansion, the current magnetic field 
is about 10~ 25 G. 



DISCUSSION AND SUMMARY 

In this letter, we discussed generation of magnetic field 
from cosmological perturbations. Separate treatment of 
protons, electrons and photons, and appropriate evalua- 
tion of collision term allowed us to obtain the evolution 
equation for magnetic field. We saw that electric current 
would be induced by collision term between electrons and 
photons. Up to the second order, there are three con- 
tributions which act as sources for magnetic field. In 
addition to well-known vorticity effect, we found a con- 
tribution from the anisotropic stress of photons. Then 
we estimated the magnitudes of two of the three con- 
tributions which can be evaluated from the first-order 
quantities. We showed the two contributions are com- 
parable at all order and the amplitude of the produced 
magnetic field is about 10 _19 G at lOMpc comoving scale 
at the decoupling. Concerning these two contributions, 
power spectrum of the produced magnetic fields and cor- 
relation with the cosmic microwave background must be 
studied. We will present them elsewhere soon [2lj . 

To evaluate the vorticity term, second-order perturba- 
tion theory is necessary. Although many authors have 
been trying this [221 1231 124I Ejf , the complete formula- 
tion including matter perturbations is still a challenging 
problem. 

K. T. and K. I. are supported by Grant-in- Aid for JSPS 
Fellows. 
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